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Abstract 

In [F. Jiang, S. Jiang, On instability and stability of three-dimensional gravity driven viscous flows 
in a bounded domain, Adv. Math., 264 (2014) 831-863], Jiang et.al. investigated the instability 
of Rayleigh-Taylor steady-state of a three-dimensional nonhomogeneous incompressible viscous 
flow driven by gravity in a bounded domain O of class C 2 . In particular, they proved the 
steady-state is nonlinearly unstable under a restrictive condition of that the derivative function 
of steady density possesses a positive lower bound. In this article, by exploiting a standard energy 
functional and more-refined analysis of error estimates in the bootstrap argument, we can show 
the nonlinear instability result without the restrictive condition. 
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1. Introduction 

The motion of a three-dimensional (3D) nonhomogeneous incompressible viscous fluid in the 
presence of a uniform gravitational field in a bounded domain O C M 3 of C 2 -class is governed by 
the following Navier-Stokes equations: 

( p t + v • Vp = 0, 

< pv t + pv ■ X7v + Vp = pAv - gpe 3 , (1.1) 

I divu = 0, 


where the unknowns p := p(t,x), v := v(t,x) and p := p(t,x ) denote the density, velocity and 
pressure of the fluid, respectively; p > 0 stands for the coefficient of shear viscosity, g > 0 for the 
gravitational constant, 63 = (0, 0,1) for the vertical unit vector, and —ge 3 for the gravitational 
force. I11 the system (11.11) the equation (11.11) 1 is the continuity equation, while (1 1.1 j) o describes 
the balance law of momentum. 

We studied the instability of the following Rayleigh-Taylor (RT) steady-state to the system 
(11.11) as in 16]: 

v(t,x) = 0 and Vp = — gpe 3 in O, (1.2) 


where 


p G C 2 (0), inf (p(rc)} > 0 and d X3 p > 0 for some xo G D. 

a;Sf2 


(1.3) 


It is easy to show that the steady density p only depends on x 3 , the third component of x. Hence 
we denote p' := d X3 p for simplicity. Moreover, we can compute out the associated steady pressure 
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p determined by p. The third condition posed on p in (11.31) means that there is a region in which 
the RT density profile has larger density with increasing X 3 (height), thus leading to the nonlinear 
RT instability as shown in Theorem 11.11 below. RT instability is well known as gravity-driven 
instability in fluids when a heavy fluid is on top of a light one. 

To investigate the RT instability of the system (II. ip around the steady-state (11.21) . we denote 
the perturbation by 

g = p — p, u = v — 0 , q = p — p, 
then, ( g,u,q ) satisfies the perturbed equations: 


g t + u ■ X7(g + p) =0, 

(g + p)u t + (g + p)u ■ Wu + Wq = pAu - gge 3 , 
divw = 0 . 


(1.4) 


(e,u )\ t =0 

= (g 0 ,u 0 ) in H 

(1.5) 

u an = 0 

for any t > 0 . 

( 1 . 6 ) 


To complete the statement of the perturbed problem, we specify the initial and boundary condi¬ 
tions: 


and 


Moreover, the initial data should satisfy the compatibility conditions Wo|<9n = 0 and divu 0 = 0. 
If we linearize the equations (11.41) around the steady-state (p, 0), then the resulting linearized 
equations read as 

gt + p'u 3 = 0, 

pu t + Vq = pAu - gge 3 , (1.7) 

divu = 0, 

where u 3 denotes the third component of u. 

Here we briefly introduce the research progress for RT instability of continuous flows, please 
refer to 12, 13, 22, 24] for incompressible and compressible stratified fluids, and 0 m 0 
for stratified MHD fluids. Instability of the linearized problem (i.e. linear instability) for an 
incompressible fluid was first introduced by Rayleigh in 1883 [23J. I 11 2003, Hwang and Guo 
[l5| proved the nonlinear RT instability of ||(p, w)IU 2 (0) i n the sense of Hadamard for a 2D 
nonhomogeneous incompressible inviscid fluid (i.e. p — 0 in the equations (11.41) 1 with boundary 
condition u-n\gQ = 0, where = {(xi,x 2 ) 6 R 2 | —l< x 2 < m} and n denotes the outer normal 


vector to dQ. Jiang et.al. 0 showed the nonlinear RT instability of 11 '^ 3 11 7^2 ^ 3 ^ for the Cauchy 
problem of (11.41) in the sense of Lipschitz structure, and further gave the nonlinear RT instability 
of ||u 3 ||l 2 ( 0 ) bi 18] in the sense of Hadamard in a unbounded horizontal period domain H. 

Recently, for a general bounded domain H, Jiang et.al. showed that the steady-state (11.21) 
to the linearized problem (ll.4l) - (ll.6l) is linear unstable (i.e., the linear solution grows in time in 
by constructing a standard energy functional for the time-independent system of (11.71) 
and exploiting a modified variational method. Based on the linear instability result, they further 
showed the nonlinear instability of the perturbed problem (jl.4[) — (jl. 6 |) by a bootstrap technique 
under the following restrictive condition (i.e., the derivative function of steady-density enjoys a 
positive lower bound): 


inf- {p\x)} > 0 . 

x£il 


( 1 . 8 ) 
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The bootstrap technique has its origins in the paper of Guo and Strauss [l0|,Lll|]. It was developed 
by Friedlander et.al. Q, and widely quoted in the nonlinear instability literature, see [J, J[@- 
[9, 21] for examples. However the Duhamel’s principle in the standard bootstrap argument can 
not be directly applied to show the nonlinear instability of the problem (ll.4jl — (11. 6 |) . see jTsf for 
the details. To circumvent this obstacle, Jiang et.al. used some specific energy error estimates 
to replace Duhamel’s principle, in which the key step is to deduce an error estimate for ( g d ,u d ) 
in L 2 (Q) (i.e the L 2 (H)-norm of difference between a nonlinear solution ( g s ,u s ) to the problem 
m-m an d a linear solution (p a , u a ) to the problem (jl.5[1 — (II .7j) J in the bootstrap technique. 
To this purpose, they introduced a new energy functional under the condition (II. 8 ft to avoid the 
integrand term 


< ((e* + p)w?)t, m? > dr, 


(1.9) 


Jo 

since the energy estimate of Gronwall-type (see (j2.3[) J does not directly offer any estimate for the 
term ((p 5 + p)w d ) T . Here < •, • > denotes the corresponding dual product between the two spaces 
and HlifiV), and represents the dual space of := {u G | divw = 0}. 

Using the new energy functional, they can get a sharp growth rate A of any linear solution (p, u ) 

in the norm “ ^/ll^llz , 2 + II' u IIl 2 (o)” • Thus, applying this property to the process of specific energy 
error estimates, they easily obtained the desired error estimate, and thus showed the nonlinear 
instability. 

This article is devoted to canceling the condition (II. 8 ft in the proof of nonlinear instability 
in [l(|. More precisely, we establish the following improved result by using a standard energy 
functional and more-refined analysis techniques to deduce the error estimate for ||(p d , M d )||L 2 ( 0 ) 
in the bootstrap argument, which will be showed in Section [3] 


Theorem 1.1. Assume that the steady density p satisfies (11.31) . Then, the steady-state (jl .2[) of 
the system, (11.41) - (11.61) is unstable in the Hadamard sense, that is, there are positive constants A, 
mo, e and So, and functions (po,uo) £ H 2 {fX) x H 2 {fX), such that for any 8 G (0, Jo) and initial 
data (po,uo) (Jpo, Jho) there is a unique strong solution ( g,u ) G C°([0, T max ), H 2 (Q) x H 2 m 
of (PHI) dLSD with a associated pressure q G G°([0, T max ), U 1 (f2)), such that 

II£ ) (T' 5 )||l 2 (T 2 ), \\(ui,u 2 ){T s )\\ L 2^, ||rt 3 (T' 5 )|| L 2 ( r2 ) > £ 

for some escape time T s := J-ln^j G (0, T max ), where T max denotes the maximal time of existence 
of the solution ( g,u ). 

By virtue of [l 6 j , the key step in the proof of Theorem 11.11 is to establish a error estimate 

|| (p d , w d ) ||i / 2 (q) < CS 3 e 3At for some constant C (1-10) 


without the restrictive condition (11.81) (i.e., Lemma 13.11) . Here we sketch the main idea in the 
proof of (ll.lOp without (11.81) . In view of the property of standard energy functional (see (I2.ip h 

A is also a sharp growth rate of any linear solution (p, u) in the norm “-y/||f?||i 2 + IMI/ryn)”’ see 

, Proposition 3.3.]. When applying the sharp growth rate of the standard energy functional 
to the process of specific energy error estimates, we need to deal with the difficulty arising from 
the term (j 1 .9 p . However, by a classical regularization method, we can show that 


16 


2 / < ((P 5 + p)u A T ) r ,u A T > dr 
Jo 

= [ (P 5 + P)\ut(t)\ 2 dx - f (p 5 (0) +p)|w d (0)| 2 da; + 
Jq J ft 



pf |ti d | 2 dxdr, 


3 




































Then, we can deduce from the error equations (see (13.91) ) that 

\Wq 5 + (t)\\ 2 L2 +2(i J IIVwJlIiadT = J gp'\u^(t)\ 2 dx + R x + R 2 {t), 

where the two higher-order terms R± and R, 2 {t) (see (j3. 14[) and (13. 15j> for their definitions) can 
be controlled by 5 3 e 3At . Using the definition of sharp growth rate, we can further estimate that 

\\VQ S + pu*{t )\\ 2 L 2 +2(i [ ||Vu£||£ a dT 

J o 

< A 2 \\VQ s + pu d m L2 + Ap\\Vu d (t )\\ 2 L 2 + C5 3 e 3At . 


Based on the estimate above, by more-refined analysis, we can further obtain the following 
Gronwall’s inequality 

^\\Vq S + PU d (t)\\ 2 L 2 + p\\Vu d (t)\\ 2 L 2 

< 2A ^|| vV + pu d (t) \\ 2 L 2 + /i J ||V-u d ||^ 2 d r^j+CS 3 e 3At . 

Since p s + p possesses a positive lower bound, thus we immediately get the desired error estimate 
(11 . 10p from the Gronwall’s inequality above and the mass equation. We mention that Jiang et.al. 
[lfi| | used another energy functional and the restrictive condition (11. 8ft to deduce the following 
Gronwall’s inequality 


d 

dt 


g d | 2 p| w d | 2 \ 

P 9 ) 


dx < 2A 


p d | 2 p|u d | 2 \ 

P 9 ) 


dx + C5 3 e 3At 


and thus got (II. 10p under (11.811 . 

Finally, we end this section by explaining the notations used throughout the rest of this 
article. For simplicity, we drop the domain f! in Sobolve spaces and the corresponding norms as 
well as in integrands over 12, for example, 


L p := L p (U), H k ■= W k ’ 2 (n), Hi := H l a (i2 ), 
In addition, we denote It ■— (0 ,T) and It ■— [0, T] for simplicity. 


in 


2. Preliminaries 

This section is devoted to introduction of two auxiliary results, which were established in [[lfij 
and will be used to prove Theorem 11.11 in next section. The first result is about the instability 
result of the linearized problem (jl.5p - fjl.7p . 

Proposition 2.1. Assume that the steady density p satisfies (11.31) . Then the steady-state (11.2 j) 
of the linearized system fll.Sj) - (ll.7j) is unstable. That is, there exists an unstable solution 

((?, u , q) := e At (-p'v 3 /A,v,p) 

to (11.51) (1 1. 71) . where ( v,p ) G H 2 x H 1 solves the following boundary problem 

f A 2 pv + AVp = ApAv + gpv 3 e 3 , 

\ div v = 0, v = 0 
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with the positive constant growth rate A defined by 


2 g J p'wjdx-Apj n \\7w\ 2 dx 


A = sup „ 

wgh- J p\w\ 2 dx 

Moreover, v satisfies v 3 fiO, vf + v 2 ^0 and 

P'v 3 £ 0 , 

where u, denotes the i-th component of v for i = 1, 2, 3. 


( 2 . 1 ) 


( 2 . 2 ) 


Remark 2.1. The linear instability was showed in 161, Theorem 1.1] except (12.2p . However, we 
can easily get (12.2p by contradiction. Suppose that p'v 3 = 0, then 


0 < A 


2 g J p'vldx — A/i f |Vu| 2 da; f |Vu| 2 da; 

/ p|h| 2 da; 


= -Ap- 


J p|h| 2 da; 


< 0 , 


which contradicts. Therefore, (12.211 holds. 


The second result is about a local existence result of a unique strong solution to the perturbed 
problem (j 1.4 p (1 1.6 p . which enjoys an energy estimate of Gronwall-type, see 16, Proposition 3.3] 
for the detailed proof. 


Proposition 2.2. Assume that the steady density p satisfies (11.31) . For any given initial data 
(q 0 ,u 0 ) G H 2 x H 2 satisfying inf xe n{/9o(®)} > 0, and the compatibility conditions u 0 \qq = 0 and 
divwo = 0, there exist a unique strong solution ( g,u ) G C°([0, T max ), H 2 x H 2 ) to the perturbed 
problem (II,4p — (11 , 6 p with a associated pressure q G C°([0, T max ), H 1 ), where T max denotes the 
maximal time of existence. Moreover, 

(1) u t e C°([0, T max ), L 2 ) and 


0 < inf {^o(^) + p} < inf {^(f, x) + p} < sup{g(t, x) + p} < sup{^ 0 (^) + p} < +00 
x ^ n x so x &n xeo 


for any t G [0, T max ). 

(2) there is a constant 5 0 G (0,1), such that if £{t) < 5 0 on some interval It C [0,T max ) ? then 
the strong solution satisfies 


£ 2 {t) + \\(u t ,A7q)(t)\\ 2 L 2 + f ||(Vu,u T , Vu r )||| 2 dr 

Jo 

< Ci (]f 0 2 + J ||(ftM)||i 2 dT]) 


(2.3) 


for any t G It, where we have defined that 


£(t) := £((g, u)(t)) = \J\\Q(t)\\ 2 L2 + ||u(t)||^ 2 , 
£ 0 := £({e,u)( 0 )) = \J\\po \\ 2 L 2 + \\u 0 \\ 2 H 2 , 
and the constant Ci > 0 only depends on p, g, p and 12. 
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3. Proof of Theorem 11.11 

Now we are in a position to prove Theorem ll.lt First, in view of Proposition 12.11 we can 
construct a linear solution 

(e\ 0 = eAt (Qo, ho) e H 2 X ( h 2 n Hi) (3.1) 

to the equations (11 ,5jl with a associated pressure q l = e At q 0 , where q 0 G H 1 , and (go,u 0 ) G 
H 2 x ( H 2 n Hi) satisfy 

||^o||L 2 ||ho3|U 2 ||(hoi,ho2)||L 2 >0, (3.2) 

£((£o;ho)) = \J 11 111^2 + ||ho||^2 = 1, 

where ho; stands for the i-th component of ho for i = 1, 2 and 3. 

Denote (qq,Uq) := 8(go,Uo), and C 2 ||(po,ho)||L 2 - Keeping in mind that the condition 
inf xe n{p(x)} > 0 and the embedding H 2 L°° , we can choose a sufficiently small h G (0,1), 
such that 

infggnipfo)} < inf (pq(x) + p{x)} for any h G (0, 8). 

Thus, by virtue of Proposition 12.21 for any 8 < 8, there exists a unique local solution ( g s ,u s ) G 
C°([0, T max ), H 2 x H 2 ) to the perturbed problem (ll.4p - fjl.6p with a associated pressure q s G 
C°([0, T max ), H 1 ), emanating from the initial data ( g 5 0 ,u 5 0 ) with £((g^,u 5 0 )) = where T max 
denotes the maximal time of existence. Moreover, 

0 < ---< mf (p (t, x) + p) (3.3) 

Z rrESZ 

and 

sup{p l 5 (t, x) +p}< sup{p 0 (x) +p}< C 3 \\g 0 \\ H 2 + ||p||loo (3.4) 

for any t G [0,T max ), where C 3 is the constant from the imbedding H 2 L°°. 

Let C\ > 0 and ho > 0 be the same constants as in Proposition 12.21 and £0 G (0,1) be a 
constant, which will be defined in (13.331) . Denote ho = min{h, ho}, for given h G (0, ho), we define 

T 8 : = —ln^|^ > 0, i.e., he AT<5 = 2e 0 , (3.5) 

A h 

T* := sup {£ G /-rmax | £((g 8 ,u 8 )(t)) < h 0 } > 0 

and 

T** := sup {t G It— | ||(p , 5 ,u 5 ) (t) \\ L2 < 2 8 C 2 e At } > 0 . 

Then T* and T** may be finite, and furthermore, 

£((g s , u s ) (T*)) = h 0 if T* < 00 , (3.6) 

|| (T**) || i2 = 28C 2 e AT ** if T** < T max . (3.7) 

Now, we denote T min := min {T 6 ,T*,T**}, then for all t G /T min , we deduce from the estimate 
(12.3j) and the definitions of T* and T** that 

£ 2 ((eV)W) + IlhWfe + [' livyill.dr 

Jo 

< Ci8 2 £ 2 ((g 0 ,u 0 )) + Cl [ \\(g 5 ,u 5 )\\ 2 L2 dr 

Jo 

< C18 2 + 4C 1 C|h 2 e 2 At /(2A) < C 4 8 2 e 2At (3.8) 
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where C 4 := C\ + 4C]C|/(2A) is independent of 5. 

Let ( g d ,u d ) = (g s ,u^) — 5(g l ,u l ). Noting that (p a ,u a ) := S(g l ,u i ) G C°([0, +oo), H 2 x H 2 ) is 
also a linear solution to (11.5j) - fll.7H with the initial data (pj), Uq) e H 2 x H 2 and with a associated 
pressure q a = Sq l G C°([0, +oo), i/ 1 ), we find that (p d , u d ) satisfies the following error equations: 


gf + p'u d = — u s ■ V g 5 , 

(g s + p)u d — pA w d + Vg d = f 5 — gg d e 3 , 
divw d = 0 

with initial and boundary conditions 

(p d (0),« d (0)) = 0, M d |an = 0 

and compatibility conditions 


(3.9) 


« d ( 0 )lan - 0 , divu d ( 0 ) - 0 , 

where we have defined that 

q A ;= q s _ g a e C°(I Tniin , H 1 ) and f s := -(/ + p)u 5 ■ Vu s - g s u*. 
Next, we shall establish an error estimate for (g d ,u d ) in L 2 -norm. 

Lemma 3.1. There is a constant C 4 , such that for all t G Ir min , 

\\(g d ,u d m\\h<C^e 3At . 


(3.10) 


Proof. Recalling that (p d ,w d ) = (g s ,u s ) — ( g s ,u s ), in view of the regularity of ( g 5 ,u s ) and 
(p a , w a ), we can deduce from (13.91 ) 9 that for a.e. t G /r min , 


^ / ie S + p)\ut\ 2 dx = 2 < ((£ 5 + p)u d ) t ,u d > - / g 5 t \u d | 2 dx 


(3.11) 


= 2 f (f t -gg t e 3 )u t dx - 2p J \Vu t | d x - j g t \u t \ dx, 
and ||\/ g s + puf ||l 2 G C°(/r min ), please refer to [2|, Remark 6 ]. Noting that 


— / p'\u d \ 2 dx = 2 p'u%v%dx, 


thus, using (13.91) 1 . we can rewrite the equality (13. lip as 

d 


dt 


[{g 5 + p)\u d \ 2 - gp'\u d \ 2 ] dx + 2 p / \Vu d \ 2 dx 


= J (2 ft + Zgu 5 ■ Vg s e 3 - g s t u d ) ■ u d dx, 

Recalling that u d (0) = 0, thus, integrating (j3.12[) in time from 0 to t, we get 

||Vg s + pi<?(()llii +2/1 f ||Vti d ||£sdr = [ gp'\uj(t)\ 2 dx +Ri + R 2 (t), 


(3.12) 


(3.13) 
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where 


and 


Ri = 


(g d + p)\ut\ 2 dx 


t =o 


Ri(t) = j J (2/r + 2gu s ■ V e s e 3 - glut) ■ ufdxdr. 


(3.14) 


(3.15) 


Next, we control the two higher-order terms R\ and i? 2 (£) above. In what follows, we denote by 
C a generic positive constant which may depend on p , g , p, A, hi and (g 0 , «o)- The symbol a < b 
means that a < Cb. 

Multiplying (j3.9[l o by u d in L 2 , we get 


(p d + p)\u d \ 2 dx = I (f - gp d e 3 + pAu d ) • u d dx. 


Exploiting (13.3ft and Cauchy’s inequality, we get 


(e 6 + P)Wt\ 2 dx < || f 5 - gp d e 3 1 || 2 + ||A 


7/ d H 2 
u 11 L 2 • 


By the definition of u a , it holds that 


lu a \\ H k = A J 5e At \\u 0 \\ H k for 0 < k, j < 2 , 


(3.16) 

(3.17) 


(3.18) 


thus, using (13.4ft . (13.8ft . Holder’s inequality and the imbedding H 2 '—>• L°°, we have 

II f 5 ~ 9 Q de 3 \\ 2 L 2 ^Il^ d |ll2 + \\(p S + P)|||cx>||m <S ||^2 + II/III2IKIII2 
<\\p d \\l 2 + S 4 e 4At , 

Noting that £ d (0) = 0, Aw d (0) = 0 and S G (0,1), chaining the estimates (I3.16K and (13.181) 
together, and taking limit for t —> 0 , we immediately obtain the following estimate for the first 
higher-order term Rp 


< 


R i =Jim j (p 5 + p)\u d (t)\ 2 dx 

112 , || A„,d/j.M |2 


Hmdle-Wlli, + ||A u »( t )||i, + 5 4 e 4A ‘) = V < 4 4 . 


(3.19) 


t —>-0 


Now we turn to estimate the most complicated higher-order term i? 2 (£)- Recalling the defi¬ 
nition of R 2 (t), we see that 

R 2 (t) = - 2 [ j [p s u a T + (, g s + p)u 5 T ■ Vu 5 + (p 5 + p)u s • Vuf] • u d dxdT 

[2gu 5 ■ A7p 5 e 3 — p s T (2 u a + u d + 2 u 6 • V« { )] • n d da;dr 
■—R2,l(t) + -R 2 , 2 (t). 

Using (13.4K . (13.8ft . (I3.17K . Holder’s inequality and the imbeddings H 2 <—>■ L°° and H 1 < -G L 4 , the 
integral term i? 2il (t) can be estimated as follows: 

ft 



R 2,l(t) 


< 


If? IU 2 ll u rrll^ 2 + 11(2 +P)IU°°ll' u 11 ir 2 11 Hr/ 1 ) || 11 n 2 dT 


< 


S 2 e 2Ar (Se AT + ||V<|| L 2 )dr 


(3.20) 


<S 3 e 3At + 


5 4 e 4AT dr 


l|Vnf||i 2 dr < S 3 e 3At . 


8 




















To estimate the second term i? 2 ,2(0; we use the mass equation (i.e. of = — fu s ■ V g 5 + p'u 3 )) and 
the formula of integration by parts to rewrite i? 2 , 2(0 as follows: 


R 


,2 ft) — J j [_(u s ■ \7 g s + p'u s 3 ) (2u* + w d + 2 u s • Vw 5 ) + 2 gu 5 ■ V^e 3 ] • u'fdxdr 
= I j \p'u 3 ( 2u * + u d + 2 u s ■ Vu 5 ) u d — 2gp s u 6 ■ Vd T u d ] dxdr 
-2 J j [ p 5 u 5 ■ V (u* + u 5 ■ Vu s ) ■ u d + p 5 u 5 ■ Vu d ■ « + w 5 ■ W 5 )] drdr 


—-^ 2 , 2,1 (f) + R2,2,2(t). 

Similarly to (13.201) . we can estimate that 


< / [IKIMIKIb + IKIb + \\A 2 h >)\\<\b + ll/IMI^MIva^Mdr 


< 


[h 3 e 3Ar (l + 5e Ar ) + (f 2 e 2 Ar ||V9 t ^|| l2 ] dr < 5 3 e 3At (l + <5e At ), 


and 


fl 2 A 2 (f) < / lk 1 ^||« 1 ^(||v<|U 2 ||<IU 2 + 


< 


+ ||<|| i 2 ||V ^|| L2 + K||^||V<|| i2 )dr 
[5 3 e 3AT (l + Se Ar ) + 5 2 e 2 Ar ||Vnf|| l2 ] dr < <f 3 e 3At (l + <5e At ). 


By the definition of £q E (0,1) in f!3.51) . 

S < de At < Se ATS < 2 for any t E /r rnin , 

Thus, summing up the estimates (I3.19p - fl3.22p . we get 

Ri + R, 2 (t) — R,\ + /? 2 ,i (t) + ^2,2,1 (£) + R2,2,2ft) < 5 3 e 3At , 
which, together with f!3 . 1 31) . yields that 

\\Vq s + put(t)\\h +2 ^ f Q tl V ^llL d r < J gp'\u%\ 2 dx + C5 3 e' 

Thanks to (12.ip . we have 

[ 9 f>'W 3 \ 2 dx <A 2 f p|u d | 2 da; + A/r f |Vu d | 2 dx 


3At 


=A 2 J (g d + p)\u d \ 2 dx + Ap J \X7u d \ 2 dx - A 2 J g d \u d \ 2 dx 
<A 2 J (g s + p)\u d \ 2 dx + Ap J \X7u d \ 2 dx + C5 3 e 3At . 
Chaining the previous two inequalities together, we obtain 

Wg S + pu?{t )\\ 2 L 2 + 2 pJ o II^H^dr 
< A 2 1| \/g 5 + pu d {t )|| 2 2 + A/i||Vu d (f )|| 2 2 + CS 3 e 3At . 


(3.21) 


(3.22) 


(3.23) 

(3.24) 


(3.25) 
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Recalling that u d G C°(lT min , H 2 ) and Vu d (0) = 0, thus, using Newton-Leibniz’s formula and 
Cauchy-Schwarz’s inequality, we hud that 

f't 


A/i|| Vw d (f)|| 2 2 = 2Ap f f d Xi up Xi u d T dxdr 

JO Jfl 1 <' , J^n 


l<i,j<3 
d || 2 


(3.26) 


< A fi ||Vu ||2 2 dr + /i / ||Vu T || L 2 dr, 

Jo Jo 

where u d T denotes the j-th component of u d . Putting (I3.25P and (I3.26P together, we have 


A 


\\Vq 5 + pu?(t)\\ 2 L2 +/i||VA 


d/ + M|2 
L 2 


< A|| yV + pu d (t )\\ 2 L 2 + 2A/i / ||Vu d ||^ 2 dr + Cd 3 e' 


3A t 


On the other hand, 
d 


di 


|| vV + pu d \\ 2 L 2 =2 / {g s + p)u d ■ u d dx + g s t \u d \ 2 dx 


<j\\VWT7hni> +M\Vq 5 + pu 


d II2 
L 2 


+ / g s t \u d \ 2 dx 


and 


g 5 t \u d \ 2 dx = — / (u d ■ Vg d + pu d 3 )\u d \ 2 dx 


/„.5M„.d|2, 


= J (2 g°u° ■ V« d — p'u\u d ) • u d da: 

< S 3 e 3M 

Putting the previous three estimates together, we get the differential inequality 


d£ 


WVq 5 + pu d (t)\\ 2 L 2 + p\\Vu d (t)\\l 2 


< 2A f|| vV + pu d 


z L 2 +p ||Vw d ||£ 2 dT +C(Pe' 


r3„3A t 


Recalling that n d = 0, thus, applying Gronwall’s inequality to (I3.28p . one obtains 


ll\/^ + 


pu d 


\ 2 L 2 +p / IIVn d ||| 2 dr < e 2At / {C5 3 e 3AT )e~ 2AT dr < 5 3 e 3At 

Jo Jo 

for all t < 7 Tmin , which, together with (I3.4p and (13 .2711 . yields that 

h d mP + IK d (*)lli» + T II Vu d ||£ 2 dr < 5 3 e 3At . 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


Finally, using the estimates (13.8p . (I3.23p and (13.3011 . we can deduce from the equations (13. 91b . 
that t 

\\g d (t)\\L 2 < [ ||£rlU 2 dr 


< 


(||w d ||^-i + || u d • Vp d || L 2)dr 


< / (5~ 2 e^ T + 5 2 e 2AT )dT<5--e^ t , 

Jo 

which, together with (13.30p . yields (13. 10p . This completes the proof of Lemma [3711 


(3.31) 


□ 
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Now, we claim that 


provided that small £ 0 is taken to be 


T d _ rri 

-L t r 


^0 


min 


5o _Cf_ m o) 
4’ 8 C 4 C 4 J’ 


(3.32) 


(3.33) 


where we have defined that mo min{||^o||L 2 7 II^osIIl 2 , ||(woi, wo 2 )||l 2 } > 0 due to (13.21) . 
Indeed, if T* = T min , then T* < oo. Moreover, from (13.51) and (13.81) we get 

S(( y g s ,u s ) (T*)) < 5e AT < 5e AJL = 2e 0 < 5 0: 


which contradicts with (13.61) . On the other hand, if T** < T min , then T** < T* < T max . Moreover, 
in view of (13. ip . (13.5j) and (I3.10p . we see that 

II(«*,«*) COHy < ||(eS,«!) (O \\ L2 + ||( e d ,« d ) (OL 

<i||(eV) (OL + v^V ^"' 2 

<SC 2 e AT " + VclPfie 3 **"/* < Se AT "(C 2 + 

<25C 2 e AT **\ 


which also contradicts with (13.7j) . Therefore, (13.321) holds. 

Since T 5 = T min , (13. lOjl holds for t = T s . Thus, we can use (I3.33jl and (13. 10p with t = T 5 to 
deduce that 

\\g S (T s )\\ L2 >||^(T 5 )|| l2 - \\g d (T s )\\ L , = 6 \\q\T s )\\ l2 - \\g d (T 5 )\\ L2 
>fc AT< llff 0 |li» - 

>2£olleolli» - Vcle o /2 > 2mo£o - Vcle if 2 > m 0 e 0 , 

Similar, we also have 

\\ui(T 6 )\\ L 2>2moeo-VcleT > rn 0 £ 0 , 

and 

||(^,^)(T 5 )|| L 2 > 2m 0 £0 - VcleT > m 0 £ o, 

where uj(T s ) denote the i-th component of u s {T 8 ) for i — 1, 2, 3. This completes the proof of 
Theorem 11.11 bv defining £ := mo£o- 
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